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ABSTRACT 


Finite-amplitude standing waves in air at ambient temperature 
contained in a rigid-walled rectangular cavity were experimentally 
investigated. The pressure waveform in the cavity was analyzed for 
harmonic content and compared to the theory of Coppens and Sanders 
as modified to include empirical losses. The waveforms for the 
nOO and OnO configurations were analyzed at strength parameters of 
about 0.5, 0.75, and 1.0, where the strength parameter is defined as 
MQ, (Y + 1)/2. M is the Mach number of the fundamental, oF is the 
quality factor of the fundamental resonance, and Y is the ratio of 
specific heats. It is found that the theoretical model accurately 


predicts the shape and approximates the amplitudes of the harmonic 


components of the finite-amplitude standing waves. 
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I. INTRODUCTION 


The purpose of this research is to make a detailed investigation 
of finite-amplitude standing waves in air at ambient temperatures in 
a rigid walled cavity, for comparison with the theory of Coppens and 
Sanders [1] modified to include observed energy losses. This in- 
ee igat on is a continuation of the experimental works of Beech [2] 
and Winn [3], with an extension to a rectangular cavity vice a long 
round tube. 

A rectangular cavity was chosen as it was a step toward the 
study of two and three dimensional waves, although this investi- 
gation concerns only one dimensional standing waves. In addition, 
it was thought that a rectangular geometry would produce over- 
tones which would be less harmonic than those for the more syn- 
metric geometry of the round tube. The dimensions of the cavity 
were chosen in an effort to space out the system resonances to 
avoid degeneracy or near-degeneracy. This was accomplished with 
some success, and thus allowed for the accurate location and anal- 


ysis of the first few overtones of each of the two modes studied. 





II. BACKGROUND AND THEORY 


A finite amplitude wave is a solution of the non-linear wave 
equation that describes the real-life environment that confronts 
large amplitude acoustic waves. This finite-amplitude wave re- 
quires second order terms, or higher, to describe the waveform. 

It has been known for many years that hydrodynamic equations 
describing traveling waves in a non-dissipative medium predict a 
change in the waveform. Kirchoff C4) and Lamb C5] were among the 
early investigators in this area. [In more recent years there have 
been many theories and methods of solution proposed to describe 
this distorted traveling waveform in a non-dissipative medium. 
Amoung these are the works of Fay Gels and Blackstock L7] who pro- 
duced an extensive theoretical treatment of low amplitude non- 
linear waves, and made use of Fourier and power series 
approximations. 

When the problem was extended to dissipative mediums, it was 
shown very early (8] that waves in absorptive media exhibit finite- 
amplitude characteristics only at high strengths, different from 
waves in a non-dissipative media. More recently, an important con- 
tribution to the understanding of finite amplitude waves ina dis- 
Sipative medium was presented by Fay Foie He, Similar to a later 
work by Blackstock Lio], formulated a Fourier series solution to 
describe the finite-amplitude waveform for a plane wave of infinite 
extent. More recent theoretical approaches to finite-amplitude 


effects have been of two basic types: an approximation solution, 





Such as that proposed by Fox and Wallace fae and a perturbation 
solution like that of Fay [9]. The approximation solution is based 
on assumed attenuation mechanisms and manner of overtone population. 
The perturbation solution makes use of perturbation methods and 
Fourier series representation of the waveform Fito. 12a A per= 
tur bation analysis which was extended to the sixth order by Keck 
and Beyer [13] is of particular importance to this investigation, 

as will be evident later. The solution obtained by Keck and Beyer 
is valid near the source, unlike previous perturbation solutions 

of finite-amplitude traveling waves in dissipative media. 

Although most investigations have been confined to traveling 
waves, there has been some investigations of standing waves. For 
example, Keller [14] extracted solutions for standing waves in a 
closed tube. In that investigation, Keller assumed a non-dissi- 
pative medium which resulted in infinite amplitude waves at re- 
sonance, yet still obtained useful results for frequencies close 
to those of resonance. Coppens and Sanders oe developed an ex- 
tension of the Keck-Beyer perturbation approach, that include wall 
losses as predicted by the Rayleigh-Kirchoff theory [15]. This 
approach was carried to the tenth order by Ruff L116) atiwhven 
point it was found that the mathematics become divergent. 

Coppens [17] later expanded the theory and used a Fourier decom- 
position technique, which both Beech [3] and Winn f2] experimentally 
Investigated. It was found that the tubes that they investigated 
did not behave in the manner Bea evenniy Rayleigh-Kirchoff loss 


mechanisms. This investigation follows up the observation that 





the tubes exhibited non-Rayleigh-Kirchoff behavior, by testing the 
theory using empirically determined losses. 

In the theory of Coppens and Sanders, a one dimensional, non- 
linear wave equation with dissipative terms which correspond to 
those encountered in the discription of plane standing waves in a 
rigid-walled cavity, is developed. The wave equation in this 


theory can be Fourier decomposed to take the form 
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and contains a dispersive term and a dissipative term, both based 


on the Rayleigh-Kirchoff theory of wall losses (15), and 


55 = 1/Q, 
: ; th : 
§ = particle displacement due to n harmonic 
ec 
n $n 


b= (1 + oy ee 


Ss = specific heat of medium at constant pressure 
eo = specific heat of medium at constant volume 
W, = frequency (angular) of fundamental resonance 

0. = iy 
1 Wy 7 Wy 


W,, and W,_ are the half power frequencies above and below the 
fundamental. 


In the description of waves in cavities, Eq. 1 can be refor- 


mulated [18] to introduce the actual cavity resonances and 





generalized to include an arbitrary, empirical absorption coeffi- 


cient for each mode. The result is 
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Wd) = resonance frequency of the n overtone 
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= nL Such that pressure antinodes exist at x = O and at x=L 


AW = nW-W 
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W = frequency at which system is being driven, near W 
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p = acoustic pressure 
Pp, = density of medium 
Co = thermodynamic speed of sound in medium 
M = mach number 
= Z 
. Boos 
Py = peak pressure of fundamental component of wave 
BC = Fourler coefficient of ae harmonic component, normalized 
such that Ry = | 
and 
th : = 
ae = phase angle of n harmonic component, where 2) = 0. 





This equation can be manipulated into a general, time-independent 
form which, when evaluated at the rigid boundaries of the cavity, 
can be further simplified to produce the following form amenable to 


computer calculation: 
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where 
Ho = 
and 
tan OQ = oe or : 
n Q,72. 


As can be seen from the above equations, it 1S important to know 
the Q(N)'s and ws of the system for the input of the theoretical 
predictions calculated by the computer. It is convenient to de- 
fine a quantity, E(n), that indicates the position of oe relative 


to the classical harmonic frequencies, nW E(n) is defined by 


le a 


E(n) = (W/ (nw =i 


i) 
The values of Q(n) and E(n) are directly determined from the 
infinitesimal amplitude analysis of the cavity and are the input 


quantities for the computer calculations. 
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IIIT. EXPERIMENTAL CONSIDERATIONS 


A. APPARATUS 

A block diagram of the experimental apparatus is shown in 
Figure 1. Finite-amplitude waves were generated in air at ambient 
temperatures in a rectangular cavity with interior dimensions of 
6.97 cm. by 20.96 cm. by 30.28 cm. The walls of the cavity were 
constructed of 0.985 in. thick aluminum plates, milled and screwed 
together to provide tight, right angle joints. These joints sealed 
by applying thin layers of silicon grease to adjoining surfaces 
before screwing them together. Acoustic waves were generated in 
the cavity by exciting the air with a plane-faced piston. The 
piston face gained access to the interior of the cavity by way of 
an orifice bored in one of the cavity walls. The piston was fitted 
with two O-rings and the piston orifice was bored such that the 
piston with O-rings would make an airtight fit. It was found that 
the piston motion was restricted by the second O-ring; therefore 
the piston was inserted into the orifice only as far as the first 
O-ring. Lubrication and further sealing of the piston-orifice 
interface was provided by the same high vacuum silicon grease used 
to seal the wall joints. The piston was attached to an M-B Elec- 
tronics Model EA 1500 Exciter. This exciter was powered by two 
M-B Electronics Model 2120MB Power Amplifiers operating in parallel. 
The maximum attainable acceleration of the piston with this equip- 
ment was about 65 g's in free air at 100 Hz. The frequency and 


power level of the piston were controlled by a General Radio 1161-A 
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Coherent Decade Frequency Synthesizer. Frequencies can be selected 
to within less than 0.01 Hz with this synthesizer. The piston mo- 
tion was monitored by measuring the output of an Endeveco Model 2215 
Accelerometer mounted within the piston. Distortion in the piston 
motion, as determined by the accelerometer, set the lower limit to 
which reliable measurements of the non-linear distortion in the 
pressure waveform could be made. Measurement of the pressure wave- 
form in the cavity was made in the far corner from the piston 
source. Refer to Figure 2. The sound pressure levels and harmonic 
components of the pressure wave in the cavity were monitored by a 
1/4 inch diameter Bruel and Kjaer type 4136 Condenser Microphone. 
This microphone was mounted such that the active monitoring element 
of the microphone was flush with the interior of the cavity. The 
microphone was fitted with an O-ring and seated in its orifice with 
the same high-vacuum silicon grease. The output signal of the 
microphone was amplified by a Bruel and Kjaer type 2801 Power 
Supply. This power supply with the microphone comprized the mi- 
crophone system. The microphone system was monitored by a Hewlett- 
Bolekard Model 400D Vacuum Tube Voltmeter. Another vacuum tube 
voltmeter of the same type was used to monitor the output of the 
accelerometer. A Hewlett-Packard Model 302A Wave Analyzer was 
used to measure harmonic distortion of both the accelerometer and 
microphone output signals. A Hewlett-Packard Model 521A Electronic 
Counter was used to monitor the filtered output of the wave ana- 


lyzer to assure that the proper frequency was being analyzed. 


Ee: 
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B. EQUIPMENT CALIBRATION 

Absolute calibration of the microphone system response was 
required to determine strength parameters. The microphone system 
included the microphone, its preamplifier, and the vacuum tube 
voltmeter connected to the output of the preamplifier. The micro- 
phone system response is defined by 

oe Vere 

where Ya is the RMS voltage from the preamplifier as indicated on 
the VIVM and P is the RMS pressure amplitude at the microphone 
face. Calibration of the microphone system was accomplished with 
the use of a Breul and Kjaer Model 4220 Pistonphone which produces 
pees + O.2 re volt/Ubar acoustic signal at 250 Hz. The micro= 
phone system sensitivity was found to be 


3 


es = (1.46 os One) 10° volt/(Nt/m*), 


which is in close agreement with 


S ee O et O sO i 1073 Volt /(Nt/m- } 


m 


and 


3 Volt/(Nt/m*) 


S 
m 


1.52 + 0.03 x 10° 
as determined by Beech [2] and Winn [3] respectively. The micro- 
phone system had a flat frequency response to within + 0.5 db 

from 80 Hz. to 20 kHz., as determined from the microphone response 
curve supplied by the manufacturer. With the microphone system 
calibrated, the wave analyzers were calibrated with respect to the 


microphone system VIVM. This was accomplished by removing the 


microphone and its preamplifier from the VTIVM, and substituting 


ie 





a voltage signal from the frequency synthesizer to both the VIVM 
and the wave analyzers simultaneously. 

The calibration of the accelerometer was accomplished by 
direct measurement of the piston face displacement at 100 Hz. and 
at various power settings with a stroboscopic light and a traveling 
microscope. The frequency of the strobe light was set very near 
100 Hz so that the beats of the piston/strobe light interaction 
could be observed. These beats, which occur at the difference 
fx equency oe ae ~ econ » allowed the piston travel to be 


easily measurable with the traveling microscope. Twelve deter- 


minations of the accelerometer sensitivity resulted in 


a G2 o0r O20) © 1074 Vole seca 
This compares with the value 
Sa = (2.35 be O03) 1074 Veneer 


determined by Beech [3]. 


©. MICROPHONE CUITPUT AND STRENGTH PARAMETER 

In order to investigate the pressure waveform in the cavity, 
as measured at the microphone, it was necessary to relate the 
strength parameter to observable quantities. “The strength para- 
meter is defined as 


See MbQ, - 


The peak pressure P, is related to the RMS pressure i by 


1 ms 


=\ 2 : 
a Pems 
It is known from the microphone sensitivity that 


We ms i: WE 
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where V is the RMS voltage reading on the microphone system VTVM. 


The strength parameter can now be reformulated as 
es Z 
sP = {2 VbQ,/(S_ pc"), 


all of which are observable or calculable quantities. The calcu- 
lation or measurement of each of these quantities introduces error 
into the determination of the strength parameter; this error is 
estimated to be about 1.0 to 1.5 percent for both the modes 


investigated. 


D. FREQUENCY PARAMETER 
The frequency parameter is a quantity which indicates the 
position of the frequency being analysed relative to the fundamental 


resonance of the system. The frequency parameter is defined by 
FP = 2(W-W, | )Q, 7. 


and can be shown to be related to the experimental frequency 


increment by 


FP = 2(W-W, )/ (W.-W _) : 


E. SYSTEM ALIGNMENT 

The lower limit to which significant measurements could be 
taken was determined by the distortion introduced by the piston 
motion. In an effort to obtain meaningful data at the higher 
harmonics and for any harmonics at low strengths, great care was 
taken to reduce the distortion introduced into the cavity. The 
distortion of the piston motion when free of the cavity was below 
the general electronic noise level of the accelerometer system and 


thus was insignificant. The bulk of the distortion came from 


ay 





misalignment of the piston and cavity orifice. The procedure used 
. to obtain maximum alignment of the piston and thus minimize the 
distortion was as follows: The output signal of the accelerometer 
system, which measured the piston motion, was displayed on the 
dual trace oscilloscope and this waveform was compared to a pure 
sine wave. The piston and exciter apparatus, which was mounted in 
an adjustable framework, was positioned relative to the cavity 
orifice until the accelerometer waveform was as clean and undis- 
torted as possible. In order to minimize any effects which may 
have arisen from vibration coupling of the cavity and exciter or 
cavity and table, the cavity was mechanically isolated from both 
the exciter and the table. This was accomplished by mounting both 
the cavity and the piston-exciter apparatus on platforms con- 
structed of alternate slabs of wood and rubber. This restricted 
the interaction of the cavity with the piston to the piston ori- 
fice only. With proper alignment, it was possible to reproduce 
conditions where the distortion in the input waveform to the 
cavity was less than 0.2% second harmonic and 0.3% for third har- 
monic, as measured by the H-P 302A wave analyzer. Distortion 
values for the higher harmonics were below the general ambient 


noise of the total apparatus, and thus were ignored. 
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IV. “PROCEDURES 


A. PRERUN PROCEDURES 

Before any meaningful data could be taken, it was important to 
determine the amount of vibration transferred through the cavity 
walls and picked up by the microphone. To make this determination, 
the top of the cavity was removed, thus relieving the pressure seal, 
and a wood plate added between the piston and microphone to isolate 
the microphone from the pressure waves induced in the air by the 
piston. As described in the System Alignment section, the cavity 
and the exciter-piston assemblies eee: isolated by a floating 
platform technique. The vibration picked up by the microphone was 
measured at frequencies of interest and produced microphone system 


5 4 


responses of 8-10 ~ to 2°10 ° volts. In all cases these values 
amounted to distortion well below 0.3%, which is the maximum 
amount of piston motion distortion. 

It was noted by Beech Eee that temperature changes produced 
significant changes in the resonant frequency of the tube he in- 
vestigated. The temperature effect on the cavity investigated 
here waS very pronounced. This presented many problems as the 
frequency of fundamental resonance was the reference for all de- 
terminations of frequency parameter. Therefore it was necessary 
to know the value of the fundamental resonance at all times. To 
minimize this shift in the resonance frequency, the entire ap- 


paratus was allowed to warm up for a period of one to four hours, 


depending on the strength parameter to be investigated, prior to 


19 





takang data. This procedure did not totally relieve the effect of 
the frequency shift so another calibration method was employed. 
Many determinations of the fundamental frequency were made during 
the taking of data, thus producing a calibration curve, see 

Figure 3, by which the frequency parameter values could be 
corrected. 

Immediately prior to the taking of each data set, the align- 
ment of the piston and cavity was checked as outlined in the system 
alignment section of Experimental Considerations. From this point 
onward, the procedure for data collection was dictated by the 
quantities being measured. There was one experimental method to 
measure the Q's and E's of the system in the linear region, and 
another to measure the harmonic content of the non-linear finite 


amplitude waves. 


B. MEASUREMENT OF Q(N) AND B(N) IN THE LINEAR REGION 

Data collection in this region required very little power, 
but required very fine control of the power setting. This control 
was obtained by inserting a voltage dividing circuit between the 
frequency synthesizer and the power amplifiers. This voltage di- 
vider reduced the available power by a factor of ten and increased 
the sensitivity of the power control by the same factor of ten. 
With this arrangement, a satisfactory fine tuning power control 
was obtained at low strengths. 

With the power set at some nominal low value, the frequency 
region about the resonant or harmonic pressure peak of interest 


was scanned until a maximum in pressure was indicated by the 
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microphone system response. The input power, determined in the 
linear, low strength, region, and kept constant during the fol- 
lowing measurements. The frequency of the driving piston was 
changed above and below the frequency of maximum pressure, until 
the pressure was -3 dB of the maximum pressure value. This deter- 
mination of the -3 dB points was made with the 302A wave analyzer, 
which measured the pressure amplitude of the cavity by way of the 
microphone system. 

It was found that the first seven overtones of the nOO mode 
and the first five overtones for the OnO mode were undisturbed 
by overlapping Q curves. The value of Q for each response curve 
was determined from the frequency difference of the -3 dB points, 
divided by the center point frequency between the two -3 dB 
points. This center point was assumed to be a good estimate of 
the peak pressure frequency, oe for the response curve in 
question. It was with these values for the system resonance fre- 


quencies that the E(n)'s were determined. 


C. MEASUREMENT OF THE HARMONIC CONTENT OF FINITE AMPLITUDS 

STANDING WAVES 

In order to obtain finite amplitude waves in the cavity, it 
was necessary to go to much higher strengths than were required 
for the measurement of the Q's and E's. With the driving fre- 
quency set to some value near the fundamental, and the peak pres- 
sure component of the fundamental held constant, the strength of 
each harmonic was measured by one of the wave analyzers. The 


other wave analyzer monitored the fundamental pressure Py which 


Ze 





was kept constant by varying the power setting of the synthesizer 
as necessary. Values of the frequency parameter were chosen such 
that the distortion curves, Figures 4 through 9, could be deter- 
mined in the areas of maximum distortion. The order in which 
measurements were taken for each distortion curve was randomized 
to avoid systematic error. Strength parameters of about 0.5, 0.75, 
and 1.0 were investigated for each mode. Strengths below 0.5 did 
not produce high enough distortion to be measured with the equip- 
ment available. Strengths above 1.0 were not accessable due to 


limitations in power. 
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Veo “RESULTS 


Tables I and II are the data for the linear, infinitesimal 
amplitude region. The information on the empirical losses is con- 
tained in the Q's and E's. These are the values used in the com- 
puter program to predict the harmonic distortion curves in Figures 
4 through 9. The theoretical curves in these figures are plotted 
along with the experimentally measured values of harmonic dis- 
tortion. The theoretical predicted values were generated for fre- 
quency parameter intervals of 0.2 and were plotted as solid curves. 
The experimentally measured values were plotted as dots where they 
occured. Data were taken, and theoretical predictions made, for 
strength parameters of 0.513, 0.782, and 1.00 for the nOO mode and 
om507, 0.780, and 1.00 for the OnO mode. 

The point of intersection of the error bars on the P,/Py theo- 


retical curve in Figure 4 corresponds to the value of EP = 2Q,E(2) 


oo 
Ne : s , : 
where W = 3 It 1S important to note that the n = 2 distortion 
curves peak when the system 1S driven at this frequency. That is, 


when the driving frequency Wis equal to the classical fundamental 
of.W,; there 1S maximum content of P,- The point where this oc- 
curs for each P,/P, curve is indicated by arrows in Figures 5 
through 9 and by crossed error bars in Figure 4. The error in the 
determination of these points comes mainly from the error in E(2) 
that appears in the above equation for FP. The error in Q) ls of 
much less significance as is indicated in Table III. The error in 


the oy curves in Figures 5 through 9 are approximately the same 


i 


as those indicated in Figure 4. 
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Table III contains the results of five separate determinations 
of the Q's and E's taken in the infinitesimal amplitude region. 
These measurements were made to provide an estimate of the error 
in the theoretical curves. In all cases the values used to gene- 
rate the curves were very close to the average values determined 
later and the regeneration of the curves was therefore unwarranted. 


Table III contains both sets of values for direct comparison. 
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Table III 


nOO Mode 
Values Used to 
Average Values Determine Curves 
3 6 
n Q(n) E(n)°10 Q(n) E(n)-10 
1 358.2 + 0.4 — ----- 357.2 |. -eee- 
2 498.5 + 4.7 6.8 + 0.1 HO Jak 720 
3 SES ts re 6 On 561.6 120 
OnO Mode 
1 Boo soe 5.3 9 eeeH- 200n Gam | ware = 
2 poe tel 256 + On 1 509.1 255 
3 M5Oneet 7.5. Sell see Oe 442.8 -11.3 
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VI. CONCLUSIONS 


It appears that the present theoretical model successfully 
predicts the major features of the harmonic content for finite- 
amplitude standing waves in a rigid-walled cavity. In all cases 
investigated, the shape and general magnitude were predicted; see 
Figures 4 through 9. 

It 1S interesting to note the position of each theoretical 
curve relative to its corresponding experimental curve. For the 
nOO mode, the theoretically predicted values are consistently 
above the measured values. The situation for the OnO mode is the 
same for P3/Py but opposite for P,/P, - A possible source of this 
difference is the large difference in the geometry encountered by 
each mode. The nOO mode is a standing wave whose wave vector is 
perpendicular to the wall with the microphone orifice, which must 
be considered a very small perturbation due to its small size. On 
the other hand, the OnO mode 1s a standing wave with its wave 
vector perpendicular to the piston orifice, which is a large per- 
turbation in the geometry. In addition, the fact that the wave 
generating surface, the piston face, 1S normal to the nOO wave 
Pector and perpendicular to the OnO vector could cause a con- 
Sestant difference in the behavior of the modes. It is certainly 
unclear to this investigator just what the specifics are that 
explain the trends in the theory/experiment relationship. [It is 
clear however, that the theory of Coppens and Sanders as modified 


to include empirical losses, does predict the essence of the 
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behavior of finite-amplitude standing waves in the rigid-walled 


cavity investigated. 
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APPENDIX A 


In the investigation by Winn geal it was found that the tubes 
analyzed did not behave in a manner consistent with Rayleigh- 
Kixchoff wall losses. [In an effort to determine if the interior 
finish of the tubes was a critical factor that determined the loss 
mechanisms, the six foot tube was reanalyzed after refinishing the 
interior to a highly polished state. These results were compared 
to ST AeonRee for the tube in a smooth but not polished condition 
as determined by Winn. 

It was found that the absorption coefficients, defined as 

WwW = WwW 
n 


_ nt 
ca > Z2C ; 
O 


were not significantly changed from the values found by Winn, and 
both sets of values obeyed the Rayleigh-Kirchoff predictions 
indicated by the solid line in Figure 10. The behavior of the 
overtones (for the polished tube) was found to be significantly 
different (from the overtone behavior in the unpolished tube), 
and was significantly different from the Rayleigh-Kirchoff pre- 


dictions as indicated by a solid line in Figure ll. 
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